Astronomy & Astrophysics manuscript no. BLCl 1 


©ESO 2012 


March 1, 2012 





Introducing the Dirac-Milne universe 

A. Benoit-Levy' and G. Chardin^ 

' UPMC-CNRS, UMR7095, Insdtut d' Astrophysique de Paris, F-75014, Paris, France, e-mail: benoitlOiap . fr 
- CSNSM, Universite Paris-Sud, CNRS-IN2P3, 91405 Orsay, France, e-mail: chardin@csnsm. in2p3 . fr 



March 1, 2012 



ABSTRACT 



(N 

o 

tin 

o^ 
(N 



o 

u 

Oh 
I 

o 



(N 
> 

in 
o 



> 

X 



The ACDM standard model, although an excellent parametrization of the present cosmological data, contains two as yet unobserved 
components, dark matter and dark energy, that constitute more than 95% of the Universe. Faced with this unsatisfactory situation, 
we study an unconventional cosmology, the Dirac-Milne universe, a matter-antimatter symmetric cosmology, in which antimatter is 
supposed to present a negative active gravitational mass. The main feature of this cosmology is the linear evolution of the scale factor 
with time, which directly solves the age and horizon problems of a matter-dominated universe. We study the concordance of this 
model to the cosmological test of type la supernovEe distance measurements and calculate the theoretical primordial abundances of 
light elements for this cosmology. We also show that the acoustic scale of the cosmic microwave background naturally emerges at the 
degree scale despite an open geometry. 
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1. Introduction 

The inflation-based ACDM cosmological model represents a 
great achievement in modern cosmology. Its predictions are in 
good agreement with a wide varieties of observational cosmo- 
logical tests, ranging from big-bang nucleosynthesis and cosmic 
microwave background anisotropics to distances measurements 
using type la supernovse. However, this success comes at an ex- 
pensive cost as two poorly-understood components are required 
in the theory to provide the concordance of the cosmological 
model to the observations. These components are dark matter, 
whose density is estimated to be around ~ 20% of the critical 
density, and dark energy, which accounts for as much as ~ 75% 
of the energy content of the Universe. 

The physical nature of dark energy remains to be elucidated 
but all the cosmological observations tend to indicate that it be- 
haves in a similar way as a vacuum energy, characterized by an 
equation of state p - -p, where pressure and density are equal in 
absolute value but opposite. Dark energy therefore appears to be 
responsible for the acceleration of the expansion of the Universe, 
which tends to be confirmed notably by type la supernovse dis- 
tance measurements. 

Faced with this uncomfortable situation where we under- 
stand less than 5% of our Universe, we study the unconventional 
cosmology of a symmetric universe i.e. containing equal quan- 
tities of matter and antimatter, where antimatter is supposed to 
present a negative active gravitational mass. The main conse- 
quence of this provoking hypothesis is that on large scales, above 
the characteristic length of the matter-antimatter distribution, the 
expansion factor evolves linearly with time, which is reminiscent 
of the Milne cosmology ( |Milne|193"3] ). 

The main motivation for considering the possibility that an- 
timatter is endowed with a negative active gravitational mass 
is provided by the study of |Carter| ( |1968[ ) of the Kerr-Newman 
geometry representing charged spinning black holes. As noted 
initially by Carter, the Kerr-Newman geometry with the mass, 
charge, and spin of an elementary particle such as an electron 



bears several of the features expected of the real correspond- 
ing particle. In particular, the Kerr-Newman "electron" has no 
horizon, has automatically the g = 2 gyromagnetic ratio, and 
has a ring structure of radius equal to half the Compton ra- 
dius of the electron. In addition, this geometry has charge and 
mass/energy reversal symmetries that strongly evoke the CP and 
T matter-antimatter symmetry (Chardin 1997 Chardin & Rax 
|1992[ |: when the non-singular interior of the ring is crossed, a 
second R'* space is found where charge and mass change sign 



(Carter|[r968 ). Therefore, starting from an electron of negative 



charge and positive mass as measured in the first R space, we 
find in the second space a "positron" of positive charge and nega- 
tive mass. The relation of the Kerr-Newman geometry to Dirac's 
equation and therefore to antimatter has been noted by some au- 
thors ( |Arcos & Pereira|2004|[Burinskii|2008l ). 

The introduction of negative masses is not a popular idea 
amongst physicists and here we do not attempt to demonstrate 
that antimatter has a negative active gravitational mass. We con- 
sider this question instead from a cosmological point of view, 
by studying the properties of this alternative cosmology and de- 
riving necessary conditions that this model must fulfill to com- 
ply with observational tests. Nevertheless, it is important to re- 
call that an analogous system exists in condensed matter, where 
electrons and holes in a semiconductor appear as quasiparticles 
with positive inertial mass, but where holes have an energy den- 
sity lower than the "vacuum" constituted by the semiconduc- 
tor in its ground stat e and antigravitate in a gravitational field 
(|Tsidirkovskii|l975|. 



The gravitational behavior of antimatter has yet to be investi- 
gated in laboratory experiments. Gravity experiments on antihy- 
drogen are under way ( [Drobychev et a l. 2007 , P erez et al.|2008| l 
and are expected to provide important insight into this behavior, 
hopefully within the next decade. 

The paper is organized as follows. In Sects. |2] and |3] we 
present some general properties of the Dirac-Milne universe. 
Section 4] is devoted to the study of primordial nucleosynthesis 
within the framework of the Dirac-Milne cosmology. In Sect. |5] 
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we Study the compliance of the model to type la supernovae dis- 
tance measurements. In Sect.|6j we investigate the geometric im- 
plications on CMB anisotropies and baryonic acoustic oscilla- 
tions. 



2. The Dirac-Milne universe 

We consider a model in which space-time is endowed with a 
Friedmann-Robertson- Walker (FRW) metric, given by the for- 
mula 



ds^ = df - aitf {dx^ + sink^dn^) , 



(1) 



where sink = sin;^',;^, sinh;t' according to the value k - 1,0,-1 
of the parameter k describing the spatial curvature, and a(t) is the 
scale factor For such a FRW metric, the time-time and space- 
space components of the Einstein tensor G^uv take the forms 



G': 



k + 2aa + 



(2) 



On scales larger than the characteristic size of the emulsion, the 
presence of equal quantities of matter with positive mass and 
antimatter with negative mass nullifies the stress-energy tensor 



Using Einstein equation G^v 
the equivalence 



%nGTijy, we therefore have 



o fl(f) oc t and k = -1. 



The metric of the Dirac-Milne universe then reads 
ds^ = dt^ - [dx^ + ^mh^xd^^) ■ 



(3) 



(4) 



As in the standard ACDM cosmology, the Dirac-Milne uni- 
verse has a distinctive geometry: while the standard ACDM 
model has a curved space-time and flat spatial sections, the 
Dirac-Milne universe has a flat space-time and negatively curved 
spatial sections. 

The Friedmann equation for the Dirac-Milne universe reads 



(5) 



where qq is the present value of the scale factor. Integrating this 
equation enables one to compute the present age of the Universe 
tu 



(6) 



where Hq is the Hubble constant. It should be emphasized that 
this is a strict equality, whereas in the standard ACDM model, 
the age of the Universe is only appro ximatively equal to 1/Hq. It 
should also be noted that the linear evolution of the scale factor 
solves the age problem of the Universe ( jChaboyer et al.|1998 1, 
which was a prime concern before the introduction of dark en- 
ergy, and does not affect the Dirac-Milne cosmology. 

Using the metric in Eq. ([T]), the particle horizon, i.e., the dis- 
tance a photon can travel since the origin, is given by the limit 



lim 

'0^0 



Jio 



dt' 

ait'Y 



(7) 



which diverges logarithmically in the case of the Dirac-Milne 
universe. This simple relation has profound implications as it 



signifies that any two given places in space were causally con- 
nected in the past. The Dirac-Milne universe therefore does not 
have an horizon problem, which removes the main motivation 
for the introduction of inflation theories. 

Therefore, the linear evolution of the scale factor and the 
property that the space-time is flat and spatial sections are open 
naturally solves two major problems in standard cosmology 
without requiring additional ingredients such as dark energy or 
inflation, and is an important motivation for studying in detail 
such a cosmology. 

Before presenting the properties of the Dirac-Milne universe 
and for the sake of clarity, we list the underlying hypothesis and 
necessary constraints that have been assumed so far: 

1. The existence of an efficient mechanism for matter and 
antimatter separation. We also assume that the Dirac-Milne 
universe is composed of separated domains of matter and 
antimatter. These universes have been abundantly studied 
in the context of inhomogeneous nucleosynthesis. Since 
this universe is matter- antimatter symmetric, this system is 
a percolation emulsion, meaning that the probability that a 
given domain extends to infinity is close to unity. It remains 
unclear, however, whether there is an efficient mechanism 
to produce such an emulsion. A fundamental parameter is 
the typical size of this emulsion. We define it as the ratio of 
the volume V to the surface S separating the two phases: 
L - VIS . At this stage of work, this characteristic size is a 
free parameter, which is assumed here to be a constant, that 
is constrained using primordial nucleosynthesis (cf Sect.|4|i. 

2. Antimatter has a negative active gravitational mass and there 
is a gravitational repulsion between matter and antimatter. 
This hypothesis is a necessity to avoid contact between 
matter and antimatter after cosmological recombination. 
Indeed, it has been demonstrated (Ramani & Puget 1976| 
[Cohen et al.|1998| l that in the case of a symmetric universe, 
where antimatter is assumed to have a positive gravitational 
mass, the global size of domains of matter and antimatter 
must be approximatively that of the observable universe. 
Otherwise, annihilations at the frontiers of matter and 
antimatter domains would generate a diffuse gamma ray 
emission that would be in contradiction with observational 
data. 

3. One of the main interests of the Dirac-Milne universe con- 
sists in providing a solution to the horizon problem without 
resorting to inflation. However, this in turn relies on the hy- 
pothesis that the contribution of radiation to the stress-energy 
tensor is negligible. If we attribute to radiation a positive con- 
tribution to the energy, the evolution of the universe will not 
differ significantly from the standard evolution for all epochs 
where the universe is radiation dominated, i.e. for redshifts 
above a few thousand. However, in the presence of equal 
quantities of particles with positive and negative masses, this 
might no longer be the case. A first step toward a full demon- 
stration that the stress-energy tensor of radiation averages to 
zero can be deduced from the papers of Hoyle & Narlikar] 



1964| , in the case studied by Hawking (1965) , where the 



particle masses are of both signs and in equal numbers. In 
this case, it can be shown that the stress-energy tensor of ra- 
diation is of order 1 jN compared to its value when all masses 
are of the same sign, and where is the number of mass car- 
riers (particles) in a Hubble horizon. 
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Temperature (MeV) 



Fig. 1. Age of the universe as a function of temperature for the 
Dirac -Milne universe (full line, red) and for the fiducial ACDM 
cosmology (dashed line, blue). 



3. Some aspects of the thermal history of the 
Dirac-Milne Universe 

Some simple properties of the Dirac-Milne universe are common 
to those of a purely linear cosmology studied in Lohiya et al. 
1998^ , [K^nghat et al.| ^\999\ [20^ , and |Sethi et al.| ( |2005| l. 



We briefly recall some of these properties. 



3.1. Time-temperature relation 

Using Eq. Q, it is straightforward to obtain the relation between 
the age of the Universe and the redshift, hence the temperature 



J_To 
Hq T' 



(8) 



where Tq is the present temperature of the Universe, as mea- 
suredby CMB experiments. To = 2.725 + 0.001 K ([Fixsen & 



[Mat her 2002). This relation between time and temperature is 
valid throughout the whole history of the universe and implies 
that the thermal history of the Dirac-Milne universe is drasti- 
cally modified from the evolution in the standard ACDM cos- 
mology. Figure ([T]l represents the age of the Universe as a func- 
tion of the temperature for the Dirac-Milne and the ACDM mod- 
els. It can be seen that, at high temperatures, the Dirac-Milne 
universe is much older than the corresponding ACDM cosmol- 
ogy. For instance, the traditional 1 MeV~l sec approximation 
for the standard model becomes 1 MeV ~ 3.3 years in the Dirac- 
Milne cosmology. As noted in Lohi ya et al.| ( |1998| l, this dif- 
ference has profound implications for big-bang nucleosynthesis 
calculations, and is discussed in section]?] 

Another temperature of interest is the temperature of the 
quark-gluon-plasma (QGP) transition. ]Omnes (197 2 ) proposed 
that matter-antimatter separation occurred around that tempera- 
ture, owing to a putative repulsive interaction between nucleons 



Fig. 2. Comparison between expansion rates and weak reac- 
tion rates. The short-dashed green line represents the neutron- 
to-proton conversion rate. The two solid lines represent the ex- 
pansion rate of the Dirac-Milne universe (red) and the standard 
cosmology (blue). The long dashed lines represent the proton- 
to-neutron conversion rates. Weak interactions decouple when 
these rates become lower than the expansion rate. 



and antinucleons. The maximum size of a domain of (anti)matter 
was controlled by the diffusion of neutrons. Aly] ( ]1974) found a 
maximum size of 7 x 10""* cm at a temperature ofT~ 330 MeV, 
which was at this epoch the estimated temperature of the QGP 
transition. This size was later found to differ from the minimum 
size a domain should have in order to ensure a production of 
primordial helium compatible with observations (Combes et al.| 
T975||Aly|197"8al l. 



The situation is rather different in the Dirac-Milne universe 
as the timescale of the QGP transition is much longer Since the 
temperature of the transition is estimated today to be around 
170 MeV ( ]Schwarz ]2003 1, it corresponds to an age of 



6 X 10^ sec in the Dirac-Mikie universe, which is a factor ~ 10"^ 
older than in the standard case. This implies that the maximum 
size to which a domain could possibly grow (assuming the ex- 
istence of an efficient separation mechanism) is five orders of 
magnitude higher. For this reason, the Dirac-Milne universe is 
far more weakly constrained by observations than the Omnes 
cosmology. 



3.2. Weak interaction decoupling 

A fundamental example of the modifications induced by a lin- 
ear scale factor can be seen in the epoch of decoupling of the 
weak interactions. This example was analyzed extensively in 
|Lohiya et al. 



1998[ l, so we only provide here the main results. 



In the standard model, weak decoupling occurs at a temperature 
r ~ 1 MeV. In the Dirac-Milne universe, this decoupling hap- 
pens at a lower temperature of around T ~ 80 keV, because of 
the slower variation and the lower value of the expansion rate. 
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This effect is illustrated in Fig.|2] Weak interactions control the 
n p equilibrium. At low temperatures, this reaction is lim- 
ited to the free neutron disintegration (green short-dashed line), 
but at temperatures higher than 80 keV, the equilibrium between 
proton and neutrons remains possible. The long-dashed line rep- 
resent the proton conversion rate as a function of the tempera- 
ture for the Dirac-Milne universe (red) and the standard cosmol- 
ogy (blue). The analytical expressions for these reaction rates 
come from Wagoner ( 1969) and Dicus et al. ( 1982 1. Weak inter- 
actions decouple when the expansion rate becomes higher than 
the p ^ n rate. The small difference in the p n rate between 
the two cosmologies is caused by a difference in the neutrino 
background temperature. As the weak interactions decouple at 
a temperature of T ~ 80 keV, neutrinos indeed also decouple 
from the photon background also at this temperature, but only 
after the annihilation of most of the electron-positron pairs. This 
implies, as noted in |Lohiya et al.| ( |1998| l, that the cosmic neu- 
trino background should have a temperature equal to that of the 
CMB. This constitutes a distinctive feature of the Dirac-Milne 
universe. 



4. Primordial nucleosynthesis 

The primordial nucleosynthesis is a key success of the standard 
model of cosmology. Theoretical predictions and observations 
are in good agreement for '^He, "'He, and D abundances, al - 
though some discre pancies exist for ^Li (see Steigman ( 2007 1, 
[Cyburt et al.| ( |2008 l for reviews of standard big-bang nucleosyn- 
thesis (SEEN)). Given the tremendous change in the timescale 
(of a factor ~ 10^ at T = 1 MeV), important modifications 



to the traditional BEN mechanism arise (Lohiya et al. 1998 
Kaplinghat et al.|1999 2000 1. Primordial nucleosynthesis in the 



Dirac-Milne universe can be described as a two-step process: 
first the thermal and homogeneous production of ^He and ^Li, 
and second, the production of D and ^He, the latter being one of 
the main novelties introduced by the matter-antimatter symme- 
try. 

4.1. Thermal and homogeneous BBN 

The age of the universe at a temperature of T ~ 80 keV, be- 
low which deuterium is able to survive its photodisintegration, 
allowing the formation of "^He nuclei, is about 40 years in the 
Dirac-Milne universe. This extremely long timescale, relative to 
the mean lifetime of the neutron, was initially considered as an 
impossible obstacle to any kind of primordial helium production 
in a linear cosmology (Kaplinghat et al. 1999). However, it has 
been pointed out ( jLohiya et al., 1998 K aplinghat et al.|[2()00 



that "^He production is possible but relies on a somewhat differ- 
ent mechanism than the standard one. 

As discussed in the previous sections, weak interactions de- 
couple in the Dirac-Milne universe at the low temperature of 
r ~ 80 keV. This implies that the neutron and proton popula- 
tions remain in thermal equilibrium until that temperature. The 
ratio of their densities is therefore regulated by the Boltzmann 
factor 



n 
P 



exp 



(9) 



where Q - 1 .29 MeV is the difference between the neutron and 
proton mass. 

For temperatures of the order of T ~ 80 keV, this number 
is a very small quantity (~ 10"^) but a few neutrons neverthe- 
less combine with ambient protons to form deuterium, which is 
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Fig. 3. ''He production rate for the Dirac-Milne (dashed red line). 
The rate drops abruptly when weak interactions decouple around 
r ~ 80 keV. As a comparison, ''He production rate for the stan- 
dard BBN is also presented (dotted blue line). 

then incorporated into helium nuclei. As the weak interactions 
are still efficient, some protons inverse-beta-decay to neutrons, 
restoring the equilibrium value. Given the timescale of 40 years, 
this very slowprocess ends up in an effective production of ''He 
nuclei. Fig. (|3 i present the production rate of "^He in both the 
Dirac-Milne universe (dashed line) and the standard cosmology 
(dotted line). We can discern the sharp decline in the production 
rate when the weak interactions decouple, leading to the nearly 
total disappearance of neutrons below T ~ 80 keV. As a sim- 
ple argument illustrating how this mechanism leads to the right 
amount of ''He, we note that the product of the ''He production 
rate with the Hubble time at production is roughly the same in 
the two models. The production rates of ''He were obtained by 
calculating the time derivative of the ''He abundances. 

To compute the primordial abundances of the different light 
elements, we solved the non-linear systems of first-order differ- 
ential equations describing the network of the nuclear reactions 



(10) 



where F, is the abundance of nuclide /, A^isot the number of nu- 
clides , A^ieac the numbers of nuclear reactions included in the 
network, f.'j the rate of reaction / + j k + I, and /^'^ the rate of 
reaction k + I ^ i + j. 

To ensure accurate results, it is necessary to include a large 
number (over a hundred) of reactions. As the time interval dur- 
ing which BBN takes place is very long, some slow reactions 
that are ineffective during SBBN must be integrated into the 
nuclear reaction network for the Dirac-Milne universe ( jLohiya 
|et al.|199"8] l. Compared to previous studies of light element pro- 
duction in linear cosmology ( [Kaplinghat et al.||2000|l, we used 
the nuclear reaction rates provided bylAngulo et al. \ 1999 1, Coc 
erni] ( [2000| , [Tang et al.| ( [2003]l,|D escouve mont et al.[^2004|i, a s 
well as the usual rates of [Wagoner, ( ,1969| )," |Dicus et al.| ( [1982| l. 



4 



A. Benoit-Levy and G. Chardin: Introducing the Dirac -Milne universe 




40 60 80 100 0.100 0.010 0.001 

rj,(, Temperature (MeV) 



Fig. 4. ^He (top) and ^Li (bottom) theoretical abundances for the 
Dirac-Milne universe as a function of the baryon density. Here 
7710 = lO'";;. Shaded areas correspond to the observational range. 



Caughlan & Fowler|([T988ll,|F ukugita & Kajino| ( [T990l l, |Rauscher 



et al.| ( |1994| l, and |Chen & Sa vage ( 1999f 



The precise theoretical predictions for the primordial abun- 
dances strongly depend - as in the standard model - on the bary- 
onic density. The final primordial abundances of "^He and ^Li as 
a function of the baryonic density, characterized by the ratio rj 
of the number of baryons to the number of photons {rj - nt/tiy), 
are presented in Fig. ([4|. The shaded areas co rrespond to obser- 
vational constraints ( [Olive & Skillman|2004l |Ryan et al.|2000| l 
on the corresponding primordial abundances. The values of rj 
such that ensure that the theoretical value of Yp is compatible 
with the observations lie in the range 

8.8 X 10"'^ < 77 < 9.6 X 10"^ (11) 

As can be seen in Fig. [4] there is, as in the standard model, 
no value of 77 that permits the compatibility of both "'He and ^Li 
with the model predictions. The lowest possible value of 77 yields 
a hthium abundance of 



H 



= 3.45 X 10" 



(12) 



Although this value is somewhat larger than the one in- 
ferred by observations, we note that it is smaller than the 
predicted value in standar d cosmology calcula tions, namely 



^Li/H= 



2008 1. The Dirac- 



(5.24+;;-^P X 10-'" ( [Cyburt et al. 
Milne universe clearly does not solve the 'Li problem, but nev- 
ertheless alleviates it. 

The value of the baryonic density inferred by the ''He ob- 
servational constraint is almost 15 times higher than the value 
usually admitted in the framework of standard cosmology. This 
high baryonic density is an important feature of the Dirac-Milne 
universe as it suppresses the need for non-baryonic dark matter 
Estimates of the dynamical mass in various structures all indi- 
cate a matter density in the Universe that is much higher than 



Fig. 5. Abundances of light elements obtained in the Milne uni- 
verse with a baryon-to-photon ratio 77 = 8.8 x 10"^. Both ''He 
and ^Li are produced at observationally compatible levels (see 
Fig. [4]i, but D and ^He are almost totally destroyed by the very 
slow thermal nucleosynthesis. 



the baryonic density deduced by the BBN in the standard cos- 
mology, which constitutes a strong motivation to postulate the 
existence of non-baryonic massive particles. In the Dirac-Milne 
universe, however, the BBN predicts a baryonic density com- 
parable to the matter density estimated by different techniques. 
There is therefore no compulsory need for non-baryonic dark 
matter in the Dirac-Milne universe. 

We note, however, that this high value of the baryonic den- 
sity aggravates the so-called missing baryon problem (F ukugita] 
2004 1. In the standard cosmology, this problem is that only ap- 
proximately half the baryons predicted by SBBN are observa- 
tionally detected. The other half is currently believed to reside in 
the warm-hot intergalactic medium (WHIM) (Bregman 2007j|. 
In the context of the Dirac-Milne universe, this WHIM would 
have to be the major source of baryons, although this possibil- 
ity remains to be investigated. We note that attempts to explain 
the dynamical behavior of galaxies consisting primarily of cold 
molecular (baryonic) gas instead of non-baryonic dark matter 
have been performed in the past (Pfenniger et al. 1994) and these 
scenarios would need to be revisited in the context of the Dirac- 
Milne model. 

Figure 5]presents the evolution of the primordial abundances 
for the light elements as a function of temperature for a bary- 
onic density of 77 = 8.8 x lO"''. It can be seen that D and ^He 
are almost totally destroyed during the stage of thermal produc- 
tion of ''He. This result was found to have severe consequences 
for linear cosmologies (Kaplinghat et al. 2000j). Although this 
is indeed the case for a regular linear cosmology without anti- 
matter domains, the presence of distinct domains of matter and 
antimatter in the Dirac-Milne cosmology provides a natural sce- 
nario for the production of D and ^He during a second stage of 
nucleosynthesis. 
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4.2. Secondary production of D and ^He 

Studies of inhomogeneous big-bang nucleosynthesis as well 
as primordial nucleosynthesis in the presence of both matter 
and antimatter domains have been conducted since the late 



70's (ICombes et al. 1975 A 


ly| 1978a Witten 1984 Alcock 


et al.||1987 Applegate et al.| 


1987 Kurki-Suonio & Sihvola 


20001 |Rehm & Jedamzik 2001 


1. It has been shown that matter- 



antimatter annihilations at the frontiers of the domains can lead 
to the production of D, ^He, and T (later decaying as ^He) mainly 
through two channels of production: nucleodisruption (p'^He and 
^He^He reactions) and photodisintegration of "*He nuclei. These 
studies provide us with the necessary material to compute the 
amount of deuterium produced by these various mechanisms. 

Our purpose in what follows is to demonstrate the possibility 
of deuterium secondary production. We therefore consider that 
the emulsion has a static behavior, in the sense that its comoving 
size is assumed to remain constant. More precise studies inves- 
tigating the dynamical behavior of the emulsion are beyond the 
scope of this first paper and will be treated in upcoming studies. 

Annihilations at the frontiers of a domain are driven by the 
diffusion of nuclei towards the frontiers. The photodisintegration 
of '^He nuclei by energetic photons resulting from electromag- 
netic cascades induced by the annihilation photons and nucle- 
odisruption are two possible processes that could produce deu- 
terium. 

The main quantity to consider is the diffusion length. It rep- 
resents the average distance over which a (anti-)nucleus can dif- 
fuse toward the frontier of the domain on a Hubble time. This 
length gives an absolute lower bound to the size of the domains, 
as any concentration of (anti)matter smaller than this diffusion 
length would be annihilated during a Hubble time. The diffusion 
length is given by |Applegate et al. ( 1987| l 



Ldiff(r) = y/6DiT)tHiT), 



(13) 



where D is the diffusion coefficient and tHiT) is the Hubble 
time at temperature T. Using the diffusion coefficients given in 
Jedamzik & Rehm ( 2001[ l and Sihvola ( 2001[ ), we computed the 
comoving diffusion length represented in Fig. |6] 

We can distinguish three regimes of diffusion. The first 
regime for T > 1 MeV is regulated by neutron diffusion. Since 
neutrons are neutral particles, their electromagnetic interactions 
with other charged particles are very weak. As the temperature 
decreases, neutrons disintegrate and diffusion is then maintained 
by protons, which being charged have a much lower diffusion co- 
efficient, which causes a dip in the diffusion length around 100 
keV. As the density decreases because of the expansion, the dif- 
fusion length gradually increases until the temperature reaches 
T ~ 50 keV. At this temperature, the density is low enough for 
the mean distance between protons to be larger than the Debye 
length of electrons. Protons therefore do not behave as free par- 
ticles but drag electrons along with them, ensuring charge neu- 
trality. These electrons are themselves subject to Thomson drag 



(Peebles 



& Rehm 



then 



1993)) and thus limit the diffusion of protons ( Jedamzik 
2001[ l. The effective diffusion coefficient of protons is 



3T 



(14) 



where crj- is the Thomson cross-section and py, the photon en- 
ergy density. 
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Fig. 6. Comoving diffusion length (solid red line). The blue 
dashed line is the comoving thermalizing length of deuterium. 
Comoving Hubble distance is also shown in green. 



4.2.1 . Annihilation rate 

In the most general case, the computation of the annihilation rate 
is a difficult task (Aly 1978b). However, with our simplified ap- 
proach based on diffusion, the estimation of the annihilation rate 
is rather straightforward. The annihilation rate is the number of 
annihilations per unit of time and surface. The quantity of matter 
(and antimatter) annihilated over a Hubble time per unit of sur- 
face is fihLdiff, so that the annihilation rate is simply this quantity 
divided by the Hubble time 



tH 



(15) 



This simple expression is found to be in good agreement with the 
one derived in Cohen et al.| fl998 ). For simplicity, and following 
previous studies of antimatter BBN (Kurki-Suonio & Sihvola 
|2000| |Rehm & Jedamzik|[200T) , we made the hypothesis that 
hydrodynamic turbulence, which could be produced by energy 
release near the domain boundary, can be neglected. 

4.2.2. Production by pliotodisintegration of "^He nuclei 

The secondary production of light elements by "^He photodisin- 
tegration can occur in many scenarios and has been extensively 
discussed in the literature. In the framework of standard co smol- 
ogy, this mechanism is known to produce D and ^He nuclei (Ellis 
let al.| 19921 iProtheroe et al.|1995| l. 

Proton-antiproton annihilations result in the production of 
neutral pions, which themselves decay to high-energy photons. 
Depending on the temperature of the background and their en- 
ergy, these photons can create e^e^ pairs on CMB -photons. 
These newly created pairs can also interact with CMB photons 
and therefore lead to the creation of electromagnetic cascades 
(Ellis et al. 1992 Protheroe et al. 1995 i. These cascades stop 
when the energy of a photon becomes lower than the pair cre- 
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Fig. 7. Quantity of D (solid green) and ^He (dashed red) nuclei 
produced by "^He nuclei photodisintegration for one pp annihi- 
lation. 



ation threshold ZSpai, = mj /Ey, where Ey is the energy of a ther- 
mal photon. Owing to the high number of thermal p hotons, the 
threshold energy for pair creation is E^ax ~ nil/22T (Elhs et al. 
1992| l. Photons with energy lower than E, 



but higher than 

Ec ~ m^/80r also undergo elastic scattering on background 
photons (Svensson & Zdziarski|1990) . The resulting spectrum 



of cascaded photons can be pa rametrized by ( Ellis et aL]|1992 
IKurki-Suonio & Sihvola|2000l l 



dny 
'dE 



[ A{E/E,) 

= A(E/E,)-' 
0, 



-^■^ E <E, 

Ec< E <E^ 
F > F 



(16) 



whereA = 3£'o£^e^/[7-(£'c/£'max)^] is a normalization constant, 
and E() is the total energy injected. 

Photodisintegration reactions ('*He(y, p)^H, '*He(y,n)-^He, 
and '*He(y,np)D) require photons with energies higher than 
the respective threshold energies 2^He(y,np)D - 26.07 MeV, 
Q4He( y.p)3H = 19.81 MeV, and e4He(y,„)-iHe = 20.58 MeV (I Cyburt| 
|et al . 2003). The existence of these threshold energies implies 
that photodisintegration is a late process, as it becomes efficient 
only when fimax becomes higher than one of the previous thresh- 
old values. This happens at a temperature Tph ~ 0.5 keV. 

High-energy photons responsible for "^He photodisintegra- 
tion mainly interact with ambient electrons by means of 
Compton scattering or with ambient nuclei by means of the 
Bethe-Heitler process (Jedamzik 2006[ ). Taking this into ac- 
count, the number of D nuclei produced per pp annihilation is 
given by 



dE. 



n„(T 4He(y,np)D(^y) d^y 



f(Ey) 



dE^ 



(17) 



where f{Ey) = npO-BniEy, 1) + riaCTBHiEy, 2) + k{Ey)neCrKN{Ey), 
(Tbh is the Bethe-Heitler process cross-section, ctkn is 
the Compton scattering cross-section in the Klein-Nishina 



4/3[ln(2£y/m,) ■ 



T979] l, and k{Ey) 



1 



regime (Rybicki & Lightman 

1/2]-' is the mean fractional energy loss by 
11995). A similar formula 



Compton scattering (Pr otheroe et al.[ 

exists for ^He. Figure |7]presents the results of this calculation. 

When photodisintegration is most effective (around T ~ 
100 eV), the mean free path of high energy photons relative to 
these processes is longer than the diffusion length, implying that 
D and ^He nuclei produced by photodisintegration will be able 
to survive and add to the overall light element production. More 
precisely, one can estimate the temperature (or equivalently the 
redshift) at which newly produced nuclei will not be annihilated 
later on. It is then necessary to assume that there exists a red- 
shift Zend where matter and antimatter cease to annihilate. The 
hypothesis of gravitational repulsion between matter and anti- 
matter leads to this gravitational decoupling but the exact mech- 
anism, which would be analogous to the separation of electrons 
and holes in a gravitational field, still needs to be established. We 
emphasize that this decoupling is a necessary condition for the 
viability of the Dirac-Milne model. A theoretical determination 
of Zend appears to be possible by analogy with the electron-hole 
system ( Ts idil' ko v skii 1 97 5 ) , but we assume here that it is a free 
parameter that will be later constrained by observations. We de- 
fine L~^°^(z) the comoving value, calculated at a redshift z, of the 
diffusion length at redshift Zcd to be 



4°ff(z) =idiif(Zend) Y 



1 +Z 
+ Zend 



(18) 



Any nuclei produced at a distance from the domain boundary 
smaller than this length will be annihilated by Zend, whereas nu- 
clei produced farther away will survive once annihilation stops. 
We can then determine the redshift z* below which the mean free 
path of high energy photons is longer than L^j'Jj!(z»), where z* is 
then the redshift below which D and ^He nuclei produced by 
photodisintegration will add up to the final abundances. 

Multiplying Eq. (17i by the annihilation rate and integrat- 
ing between the times corresponding to z* and Zend, we get N^, 
the number of D nuclei produced by photodisintegration of "^He 
nuclei 



ND{zmz)Siz) 



dz 



(l+z)2 



(19) 



We can get rid of the surface term S (z) by dividing eq. ( [T9| by 
the total number of baryons in the volume V, ni,V, which yields 
the final D abundance 



ph 



(20) 



This expression depends only on two variables, Zend^ the redshift 
of gravitational decoupling and Likev, the comoving emulsion 
size at 1 keV. The values of the final D and ^He abundances as a 
function of these two parameters are presented on Fig. |8] 

Electromagnetic cascades initiated by annihilation photons 
cause the injection of non-thermal energy into the background 
radiation. In the standard cosmology, these energy injections are 
tightly constrained, most notably by the FIRAS measurements 
( |Fixsen et al.|1 996 ). However, in the Dirac-Milne universe, ow- 
ing to the low value of the expansion rate, the radiative processes 



that can thermalize virtually any energy injection ( Danese & de 
|Zotti|1977|rHu & Silk|19 9T) decouple at a lower redshift than in 
the standard case. The detailed calculation of the thermalization 
of CMB distortions will be presented elsewhere, but the exclu- 
sion contour in the (zend - ^ikev) plane is shown in Fig. [8] 
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10" 



Fig. 8. D and ^He final abundance as a function of the size of the 
emulsion and the redshift of gravitational decoupling. Hatched 
region is forbidden because of CMB distortions constraints. 



To ensure that the fi nal D abundance ag rees with observa- 
tions (D/H ~ 3 X 10"^, [Pettini et aL]|2008| l and imposing the 
gravitational decoupling as late as possible, we find that the co- 
moving size of the domains at 1 keV has to be around 



^IkeV 



10'^ 



cm, 



(21) 



which corresponds to a comoving size of 7 kpc today. This es- 
timate of the size of the emulsion may seem small compared, 
for instance, to the typical size of a galaxy, implying that there 
should be numerous matter and antimatter domains inside a sin- 
gle galaxy. However, we recall that this size is only constrained 
when gravitational decoupling occurs and that domains, which 
are in a non-linear mode of evolution immediately after photon 
matter-radiation decoupling, will have a dynamical behavior af- 
ter recombination that effectively increases this size (Dubinski 
[etll . 1993 Piran 1997). 

As suggested by Fig. |7] the secondary production of light 
elements by photodisintegration of "^He leads to an overproduc- 
tion of ^He relative to D. It is indeed usually assumed that D can 
only be destroyed by stellar processes. On the other hand, ^He 
can be either produced or destroyed, but the ratio D/^He can 
only decrease ( |Sigl et aL]|I995[ ). Since this ratio is observed in 
our Galaxy to be D/^He ~ I but predicted to be D/^He ~ 0.1 in 
the Dirac-Milne universe, the overproduction of ^He is a priori 
a strong constraint on the Dirac-Milne model. However, we note 
(see |Jedamzik| ( |2002| l for a review) that D enhancement could 
occur at high redshift diminishing the efficacy of the ^He con- 
straint. 



4.2.3. Production by nucleodisruption 

Nucleodisruption has been found to be a significant producer of 
D and ^He nuclei ( Sihvola 200 1 ) within the standard evolution 
of the scale factor. However, with the slow evolution of the ex- 
pansion rate in the Dirac-Milne universe and the hypothesis we 



made about the spatial repartition of the matter and antimatter 
domains, this situation changes. 

Nuclei produced by nucleodisruption possess a kinetic en- 
ergy ranging from a few MeV for nuclei to a few tens of MeV 
for nucleons ( jBalestra et al.|[T988 ). These newly produced nu- 
clei thermalize by Coulomb scattering on ambient protons and 
electrons. The thermalization length for D nuclei produced with 
an energy Eq = 10 MeV is presented in Fig. |6] (blue dotted 
line). This distance is always much smaller than the diffusion 
length, implying that any D nucleus produced by nucleodisrup- 
tion will finally return towards the annihilation zone and be de- 
stroyed there. A possible way to produce a higher fraction of 
deuterium by nucleodisruption would be to consider small do- 
mains of (anti)matter within a larger domain of antimatter (mat- 
ter). This situation occurs continually in an emulsion, which suf- 
fers a redistribution of "domains" when bridges in the emul- 
sion disappear by annihilation. If the dimension of the larger 
domain is larger than the diffusion length, then an important 
fraction of the D and ''He produced by nucleodisruption could 
survive. However, precise calculations of this production require 
the knowledge of the statistical properties of the spatial distri- 
bution of domains, which strongly depends of the separation 
mechanism. This point should be investigated in future studies 
of nucleosynthesis in the Dirac-Milne universe. 



5. Type la supernovae 



In 1998 dRiess et al.|199 8"Perlm utter et al.| 1999^, distance mea- 
surements for type la supernovse (SNe la) revealed that these 
objects are dimmer than expected if our Universe was correctly 
described by a decelerating Einstein-de Sitter model. The intro- 
duction of a cosmological constant A in the field equations of 
general relativity, which is apt to produce an accelerating expan- 
sion, provided an impressive fit to the observational data. Today, 
SNe la are one of the most important cosmological tests and are 
considered as prime evidence of an acceleration of the expan- 
sion. We recall however that the strong evidence of a recent tran- 
sition between a decelerating phase and an accelerating phase of 
expansion heavily relies on the prior hypothesis of spatial flat- 
ness. Without this hypothesis, the evidence is less clear ( Seikel 
|& Schwarz|2008| l. 

The Dirac-Milne universe has neither acceleration nor decel- 
eration and is therefore equivalent to an open empty universe. In 
terms of the usual cosmological parameters, this universe corre- 
sponds to the combination 



(Qm = 0, Qa = 0) . 



(22) 



In this context, the luminosity distance in the Dirac-Milne 
universe follows the simple expression 



diiz) 



Ho 



(1 -Hz)sinh[ln(l -Hz)]. 



(23) 



It is usually claimed that the empty universe, hence the 
Dirac-Milne universe, is strongly inconsistent with the SNe la 
observations. We wish here to elaborate on this statement using 
the data of the first release of the SNLS collaboration (Astier| 
|et al. 2006). The SNLS data consist of two distinct datasets. The 
high-redshift sample, from the SNLS, comprises 71 SNe la with 
redshifts between 0.2 < z < 1.01. The second sample is a low- 
redshift set, consisting of 44 SNe la taken from the literature 
with redshifts z < 0.15. These data come from different experi- 
ments and are therefore possibly subject to different sources of 
systematic errors. 
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Following the definition given in |Astier et al.|p006) l, the dis- 
tance modulus is 



Mb 



M + a{s- I) -/3c, 



(24) 



where M is the absolute magnitude of SNe la, a and /3 are global 
parameters that link the stretch s and the color c to the distance 
modulus, and is the apparent magnitude of the supernova. It 
should be emphasized that, in contrast to the ACDM cosmology, 
there is no cosmological parameter dependence in the Dirac- 
Milne luminosity distance. The only degrees of freedom are the 
nuisance parameters, M, a, and /3. Following the procedure de- 
scribed in |Astier et al. (20061, we minimize the expression 



(//ji-51ogio(QfJ10pc))^ 



(25) 



Here, cr(pg) takes into account measurement errors in the appar- 
ent magnitude m*g, stretch, and color parameters derived by the 
analysis of light curves ( |Guy et aL]|2005| l, cri„t is the so-called 
"intrinsic" dispersion, which is a parameter introduced to ac- 
count for SNe la being astrophysical objects that naturally have 
some intrinsic dispersion in their absolute magnitude. However, 
the value of this parameter is unknown, and in the fitting pro- 
cedure, (Tint is adjusted to ensure that the reduced chi-squared is 
unity. 

5.1. Analysis with only the high-z sample 

We first performed our analysis on the high-z sample without 
including any low-z SN la. Without this low-z anchoring, the 
analysis does not permit us to discriminate between the ACDM 
and the Dirac-Milne universes. In this respect, we note that the 
three-year analysis of SNLS using their data alone (Guy et al. 
|2010| l is consistent to a better than 68% CL with the Dirac-Milne 
universe, while the Einstein-de Sitter (EdS) model is clearly ex- 
cluded. The evidence of an expansion acceleration therefore re- 
lies on a comparison between low-z and high-z SNe la. We also 
present the results for the EdS model. 



Table 1. Values of the different parameters for the fit using only 
the 71 high-z SNe la of the SNLS. Here, 0-^^ has been fixed to 
o'int = 0. The number in parenthesis indicates the number of 
degrees of freedom. 



Parameter 


Dirac-Milne 


ACDM 


EdS 






0.289 ± 0.033 


1 


M 


-19.24 


-19.36 


-19.04 


0"mt 











total 


553.64 


558.6 


724.38 


xVdof 


8.14(68) 


8.33 (67) 


10.65(68) 



The results of the analysis of the 71 SNe la are presented in 
Tables [T] and 2] In Table [T] the intrinsic dispersion is fixed to a 
null value. The analysis is therefore performed using only the 
measurement errors, thereby giving a stronger weight to SNe la 
with redshifts 0.2 < z < 0.4, which have smaller errors. The total 
and reduced x^ of the EdS model are much larger than those of 
the ACDM and Dirac-Milne models, which are in turn similar 
In Table |2] the intrinsic dispersion is determined by assuming 
that the reduced chi-squared is unity. The value of the intrinsic 
dispersion required to renormalize the reduced of the EdS 
model to unity is here again much larger than those required 



Table 2. Values of the different parameters for the fit using only 
the 71 high-z SNe la of the SNLS. Here o-jnt has been adjusted 
such that;^'^/d.o.f. = 1. The number in parenthesis indicates the 
number of degrees of freedom. 



Parameter 


Dirac-Milne 


ACDM 


EdS 


Qui 




0.25 ± 0.08 


1 


M 


-19.18 


-19.33 


-18.87 


0"int 


0.1172 


0.1173 


0.155 


total 


68.01 


67.01 


67.97 


X^ldoi 


1.00 (68) 


1.00(67) 


1.0 (68) 



for the other two models. The conclusion of this first analysis is 
twofold. First, we confirm that the decelerating Einstein-de Sitter 
model is extremely unlikely, which should come as no surprise. 
Second, the Dirac-Milne and the flat ACDM models are almost 
identical, Dirac-Milne being in even closer agreement with the 
data than the flat ACDM model. Our analysis stresses that all 
SNe la analysis depend strongly on the use of low-z data to an- 
chor the Hubble diagram. 

5.2. Analysis on the full data set 

Similarly, we proceeded with our analysis of the full data sam- 
ple used by SNLS in its one-year analysis, i.e. using an het- 
erogeneous sample of low-z SNe la. Our results are given in 
Tables|3]and|4] In this case, as expected, the flat ACDM pro- 
vides a closer fit to the data than the Dirac-Milne universe. 



Table 3. Values of the different parameters for the fit using the 
complete data set of 115 SNe la. Here, (Tint has been fixed to 
(Tint = 0. The number in parenthesis indicates the number of 
degrees of freedom. 



Parameter 


Dirac-Milne 


ACDM 


EdS 






0.258 ±0.019 


1 


M 


-19.258 


-19.38 


-19.01 


0"i„t 











X^ total 


877.93 


809.75 


1507.83 


;r'/dof 


7.83(112) 


7.295 (111) 


13.46(112) 



Table 4. Values of the different parameters for the fit using the 
complete data set of 115 SNe la. Here (Tint has been adjusted 
such that;if^/d.o.f. = l.The number in parenthesis indicates the 
number of degrees of freedom. 



Parameter 


Dirac-Milne 


ACDM 


EdS 






0.250 ±0.036 


1 


M 


-19.217 


-19.331 


-18.98 




0.1432 


0.1289 


0.258 


x'- total 


112.05 


110.96 


111.95 


;rVdof 


1.00(112) 


0.99(111) 


0.99(112) 



In Fig.|9j we present the residuals of the Hubble diagram for 
the Dirac-Milne, the flat ACDM, and Einstein-de Sitter models. 
The left-hand column represents the residuals when the value 
of the intrinsic dispersion was adjusted to normalize the x^ to 
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Fig. 9. Residuals of the Hubble diagram for the Dirac -Milne ((a) and (b)), flat ACDM ((c) and (d)), and Einstein de-Sitter ((e) and 
(f)) models. The left column represents the residuals obtained by the minimization of the;i^^ defined by eq. (25 i, in which the values 
of the intrinsic dispersion is adjusted so that;^'-/d.o.f. = 1. In the right column, this intrinsic dispersion is fixed to 0. 



1 per degree of freedom. In the right-hand column, we present 
the residuals obtained when the intrinsic dispersion parameter 
is fixed to the value o-jnt = 0. Setting cri„, to enables one to 
consider the "real" measurement errors. It appears that for SNe 
with redshifts in the range 0.2 < z < 0.8 these errors are smaller 
than the intrinsic dispersion. The use of such an ad-hoc param- 
eter in the analysis may degrade the quality of the data in the 
redshift interval 0.2 - 0.4, where the quality of the observations 
is highest. We note that such an analysis with statistical errors 
only was previously performed by the SNLS three-year analysis 
( |Guy et al.|20T0l ). 

Panels e) and f) present the residuals for the Einstein-de 
Sitter model. Even though the reduced has been constrained 
to unity in panel e), the characteristic slope in the residuals illus- 
trates the non-conformity of the Einstein-de Sitter model to the 
data. The difference between the Dirac-Milne and the ACDM 
models is however almost negligible, and can hardly be seen by 
simply examining the residuals. 

We note that it is only the use of nearby SNe la that en- 
ables us to distinguish between the Dirac-Milne and the ACDM 
models, the latter then being the most likely model at a confi- 
dence level of more than 3cr as previously announced (Astier 
let al.|2006||Kowalski et al.|2008| . 

It appears that sources of previously unaccounted for sys- 
tematic errors ( Kelly et al.|2010| l are present in the nearby sam- 
ple data set. To investigate this, we determined the constant off"- 
set to the apparent magnitude of nearby SNe la required for the 
chi-squared for the Dirac-Milne and the flat ACDM to become 
equal. We found that an off'set of dm^ = 0.06 mag is sufficient 



to ensure that the two models are equally probable. This value 
should be compared to the budget of systematics errors esti- 
mated in recent studies: AM = 0.04 mag ( Kowalski et al.|2008[ l. 
Therefore, a relatively mild systematic error of 1 .5 cr for nearby 
SNe la would lead us to favor the Dirac-Milne universe over the 
conventional ACDM cosmology in the SNe la analysis. 

6. Other tests 

A major result of CMB experiments has been the precise mea- 
surement of the position of the first acoustic peak on the degree 
scale, which seems to imply that the spatial curvature is nearly 
zero (Komatsu et al."2010). In the open spatial geometry of the 
Dirac-Milne universe, this position would naively be expected at 
a much smaller angle. The ratio of the angular distances in the 
two models taken at redshift z ~ 1 100, which corresponds to the 
surface of last scattering surface, is 



^Milne(^) ,= 1100 



jACDM 

"a 



(z) 



« 169. 



(26) 



The value of this ratio implies that an astrophysical object at 
redshift z - 11 00 is seen under an angle 169 times smaller in the 
Dirac-Milne universe than in the ACDM cosmology. 

The angular position of the first peak is defined by the angle 
under which the sound horizon is seen at recombination 



dA(z*y 



(27) 
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where XsiZt) is the sound horizon, dAiz*) the angular distance, 
and Zt the redshift of the last scattering surface. It is of interest to 
consider the equivalent multipole ~ ^/S- The sound horizon 
is defined as the distance that acoustic waves can travel in the 
primordial plasma. Taking into account the universe expansion, 
this distance reads 



Xs 



Jo 



dt' 



(28) 



where the speed of sound Cs - c/ V3(l + R), R being a correc- 
tive factor caused by the presence of baryons ( jHu & Sugiyama 
1995| l. Its value can be related to the baryon to photon ratio rj by 
R^l.lx 10'V(1 

The definition of the lower bounds of the integral requires 
some care as we have shown previously that this integral di- 
verges near the initial singularity. However, the mechanisms of 
sound generation in the Dirac-Milne and ACDM universes dif- 
fer radically. In contrast to the generation of sound waves in 
the standard model where inhomogeneities are produced at the 
epoch of inflation, it seems reasonable to consider that sound 
waves in the Dirac-Milne universe are probably produced by an- 
nihilation at the matter-antimatter frontiers. It is therefore natu- 
ral to consider for this value the epoch of QGP transition around 
T ~ 170 MeV. In the absence of inflation, the universe is indeed 
very homogeneous before this event, and QGP transition seems 
naturally to be the scale of interest in the Dirac-Milne cosmol- 
ogy. There is at present a rather general consensus that the QCD 
transition is of neither first nor second order at null chemical po- 
tential, but is probably an analytic crossover (see e.g. |Aoki et al.| 
p006l l, [Baza"vov et al.l (20Q9 ), E ndrodi et al.1 ( |20n] i, and refer- 
ences therein). Under this hypothesis, it would be extremely dif- 
ficult, if not impossible, to understand how significant amounts 
of matter and antimatter could have survived annihilation in the 
very slow evolution of the Dirac-Milne universe. However, this 
present belief is based on a QCD calculation, whereas the sit- 
uation of the primordial universe is more complicated, includ- 
ing additional light leptons (electrons, positrons, neutrinos, and 
antineutrinos). In addition, present calculations use light quark 
masses for the u and d quarks that are significantly higher than 
the actual masses of these quarks. Therefore, the observation 
of large matter-antimatter domains would be an extremely use- 
ful indication that, in contrast to present expectations, there is a 
sharp transition allowing the survival of significant regions of an- 
timatter at the QCD transition. It should also be noted that some 
authors are clearly considering the possibility that baryogenesis 
occurs at the QCD transition or at a (9(100 MeV) temperature 
(see for example |Dolgov ( 1992| l for a review). 

Acoustic waves then propagate in the plasma as long as mat- 
ter and antimatter are in contact, i.e. until the gravitational de- 
coupling, estimated in the previous section at Zen^i ~ 3 x lO'*. 
With these values, the comoving sound horizon is found to be 

- 42 Gpc. (29) 

The expression of the angular position of the first acoustic peak 
then follows 



(30) 



It can be shown that the redshift of the last scattering surface 
in the Dirac-Milne universe is a few percent lower than in the 
standard cosmology, here again due to the late decoupling of 



the radiative processes leading to recombination. We found that 
~ 1040. 

Calculating the multipole of the acoustic scale using expres- 
sion ( [30| , we obtain fa ~ 160. The standard value of this quantity 
is {„ ~ 300 ( jSpergel et al.|2003 1. Instead of a discrepancy of a 
factor » 169 , there is an almost exact compensation between the 
larger geometrical term, induced by the open geometry of the 
Dirac-Milne universe, and the larger sound horizon, caused by 
the slow evolution of the expansion rate before recombination. 
Taking into account the numerous approximations in the model, 
this remarkable coincidence is quite unexpected and represents a 
fascinating motivation to study the Dirac-Milne universe in more 
detail. 

However, the sound horizon scale is also imprinted in the 
large-scale structure power spectrum under the form of small 
oscillations (Eisenstein & Hu 1998 ) called baryonic acoustic os- 
cillations (BAOs). These oscillations are expected, at least in a 
first approximation, on the same scale as the sound horizon. As 
discussed above, the sound horizon is much larger in the Dirac- 
Milne universe than in the standard cosmology. Admitting that 
there are BAOs in the Dirac-Milne universe, they should be ex- 
pected on a scale much larger than that of the standard cosmol- 
ogy. The claimed detection of these BAOs on the expected scale 
(within the standard cosmology) of ~ 150 Mpc ( |Eisenstein et al. 
2005 ), presently detected at the ~ 3cr level, if confirmed by on- 
going experiments, would therefore provide a strong constraint 
on the Dirac-Milne universe. 



7. Conclusion 

Since the standard ACDM model is in good agreement with 
observations but rather poorly theoretically motivated, we have 
studied here an alternative cosmological model, the Dirac-Milne 
universe. Inspired by the work of Dirac, Kerr, and Carter, this 
model restores the symmetry between matter and antimatter. 
Relying on the symmetries of the Kerr-Newman solutions in 
general relativity, it makes the hypothesis that particles and an- 
tiparticles behave similarly to quasiparticles such as electrons 
and holes in a semiconductor, and that antimatter has a negative 
active gravitational mass. A fundamental characteristic of this 
Universe is the linear evolution of its scale factor, which solves 
in an elegant way both the problems of the horizon and the age 
of the universe. 

For primordial nucleosynthesis, we have found that the 
Dirac-Milne universe is able to produce "^He at an adequate level, 
while producing ^Li nuclei in proportions admittedly a factor 
three higher than the observed values but with a smaller dis- 
agreement between observations and predictions than the stan- 
dard cosmology. We have also shown that surface annihilations 
at the frontiers of matter and antimatter naturally lead to the pro- 
duction of D nuclei, in amounts that are proportional to the in- 
verse of the characteristic size of the emulsion. The main focus 
of this study has been the production of D, although the assump- 
tion that the emulsion has a fixed characteristic size is clearly an 
approximation. Relaxing this assumption will allow the possibil- 
ity of the total annihilation of small patches of antimatter inside 
larger regions of matter. This might lead to a net production of D 
by nucleodisruption, which has a higher D/^He production ratio 
(Balestra et al. 1988[ l, hence alleviate the constraints made by the 
overproduction of ^He nuclei . 

The Dirac-Milne universe does not undergo an accelerated 
expansion and therefore seems to conflict with the usual inter- 
pretation of type la distance measurements. We have shown that 
the Dirac-Milne universe can nevertheless be reconciled with 
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type la supernovae observations if we take into account a rea- 
sonable size of systematic errors in the low-z supernovae subset. 
Finally, and perhaps the most surprising result of this unusual 
cosmology is that the acoustic scale naturally emerges at the de- 
gree scale, despite the open geometry. 

These first results are encouraging but some issues have not 
been addressed in this study. Amongst these is how stracture 
formation occurs in the presence of separate domains of positive 
and negative mass. The use of numerical simulations will most 
probably be a necessity as the usual linear approximation does 
not hold. Immediately after matter-radiation decoupling, density 
contrast is indeed on the order of unity for any distribution of 
matter with positive mass and antimatter with negative mass. 
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